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ABSTRACT 

Most of the observed extrasolar planets are found on tight and often eccentric orbits. 
The high eccentricities are not easily explained by planet-formation models, which 
predict that planets should be on rather circular orbits. Here we explore whether fly- 
bys involving planetary systems with properties similar to those of the gas giants in 
the solar system, can produce planets with properties similar to the observed planets. 
Using numerical simulations, we show that fly-bys can cause the immediate ejection of 
planets, and sometimes also lead to the capture of one or more planets by the intruder. 
More common, however, is that fly-bys only perturb the orbits of planets, sometimes 
leaving the system in an unstable state. Over time-scales of a few million to several 
hundred million years after the fly-by, this perturbation can trigger planet-planet 
scatterings, leading to the ejection of one or more planets. For example, in the case 
of the four gas giants of the solar system, the fraction of systems from which at least 
one planet is ejected more than doubles in 10 8 years after the fly-by. The remaining 
planets are often left on more eccentric orbits, similar to the eccentricities of the 
observed extrasolar planets. We combine our results of how fly-bys affect solar-system- 
like planetary systems, with the rate at which encounters in young stellar clusters 
occur. For example, we measure the effects of fly-bys on the four gas giants in the solar 
system. We find, that for such systems, between 5 and 15 per cent suffer ejections of 
planets in 10 8 years after fly-bys in typical open clusters. Thus, encounters in young 
stellar clusters can significantly alter the properties of any planets orbiting stars in 
clusters. As a large fraction of stars which populate the solar neighbourhood form 
in stellar clusters, encounters can significantly affect the properties of the observed 
extrasolar planets. 
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1 INTRODUCTION 

More than 470 planets have been discovered since the first 
disco very of an extrasolar planet around a main-sequence 
star (|Mavor fc Quelojl 19951 ). The properties of most of these 
planets are very different from our solar system, with gas- 
giant planets on often much tighter and more eccentric orbits 
than seen in the solar system (http://www.exoplanet.eu). It 
is estimated that around 10 per cent of solar mass stars host 
gas -giant planets on orbit s with semi-major axes less than 5 
au |Cumming et al.l 2008). In particular, the large spread in 
eccentricities poses a problem for planet formation, which 
is believe d to occur in discs around yo ung stars (see, for 
example. IPollack et al.lll996l : lBosslll997l ). This process will 
tend to leave planets on almost circular orbits. 
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It is today commonly believed that the observed eccen- 
tricities are caused by so-called p lanet-planet scattering, first 
discussed in iRasio fc Fori (|i 996). Several models have been 
suggested to explain why planetary s ystems undergo such 
scattering. iMoorhead fc Adams! (|2008l ) suggest that systems 
containing two giant planets, which are migrating in a pro- 
toplanetary disc, can perturb each other's orbits, exciting 
eccentricities. Such systems are unstable and often result in 
the removal of one planet due to ejection, collision or accre- 
tion onto the host star. The remaining planet is left on an 
eccentric orbit after the gas disc has dispersed. 

It is however also possible that planetary systems 
evolve past the disc phase without undergoing scatter- 
ing. If the planets are too tightly packed, such systems 
may become unstab le on time-scales of a few to several 
10 m i llion years (e.g. Juric fc Tremainel l2008l:lFord fc Rasio 
l2008l : IChatteriee et all 120081 : iRavmond et"al]|2009l ). While 
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planet-planet scattering in these types of systems can ex- 
plain the observed eccentricities, migration is needed to 
explain the observed semi-major axes (see, for example, 
iMatsumura et al.ll201ol ). 

A somewhat different mechanism to explain the on- 
set of planet-planet scattering in planetary systems is that 
they are triggered by encounters in the birth environment 
of many stars; stellar clusters. In such clusters, the high 
number density of stars causes encounters between stars 
to be common. So-called exchange encounters can leave 
planet-hosting stars, which formed single, in incline d bina- 
ries w here the so-called Kozai mechanism can operate (|Kozail 
1962). This can increase planetary eccentricities, particu- 
larly the eccentricity of the outermost planet, potentially 
triggering planet-planet scattering in multi-plan et systems 
l|Malmberg et al.ll2007j ; lMalmberg fe Daviedl2009l ) . However, 
depending on the binary fraction in stellar clusters, fly-bys 
may be much more common than exchange encounters. Es- 
sentially, in clusters with a high binary fraction, exchange 
encounters are likely to dominate, while in clusters with a 
low binary fraction, fly-bys are likely play the most impor- 
tant role. Fly-bys can, if they occur during the protoplan- 
etary disc stage, cha nge the frequency of planet formation 
l|Forgan fc Rice 200 9]) or the orbital and ph ysical properties 
of planets formed (jFragner fc Nelson! [2009 1. 

Fly-bys occurring after the stage of giant planet 
formation can excite eccentr icities in planetary systems 
( Zaka mska fc Tremaind \2004 ) and trigger planet-planet 
scatt ering in relatively packed systems |Spurzem et al.l 
120091) . 

Young stellar clusters also contain many young massive 
stars, who se strong radia t ion can damage p rotoplanetary 
discs (e.g. lArmita gej l2000l : [Ad ams et all 2006). Our under- 
standing of the birt h environmen t of the solar system was 
recently reviewed bv lAdamj (|2010h . For example, the orbital 
elements of Sedna, along with the chemical composition of 
objects in the solar system, suggest that the Sun formed in 
a cluster with roug hly 10 3 -10 4 stars. 

In this paper we measure the effects of fly-bys on plan- 
etary systems resembling the solar-system, i.e. systems with 
gas giants on initially well-spaced and circular orbits. The 
effects of such fly-bys on both individual planetary systems, 
as well as on the population of planetary systems as a whole, 
can be significant. Below we summarise the most important 
findings from our experiments. 

(i) Close fly-bys can lead to the immediate eje ction of one 
or more planets from t he host system (see also Hills! 19841; 
Hurley fc Sharal |2002| ; iFregeau et all 120061 ; ISpurzem etakl 



20091 ). Such planets may either be left unbound (ionisation) , 
or be captured by the intruder star (exchange). We measure 
what fraction of fly-bys that causes ionisation and exchange 
to occur in Section 3, and in Table [5] we present the rates 
at which such encounters occur in stellar clusters with prop- 
erties representative of the young cluster population in the 
solar neighbourhood. For example, in a cluster with initially 
700 stars and half-mass radius of 0.38 pc, 3 per cent of stars 
hosting a planetary system similar to the four gas giants of 
the solar system will suffer a fly-by leading to the immediate 
removal of one or more planets from the host star. 

(ii) Wider fly-bys, which do not immediately remove plan- 
ets from the host star, may still cause the ejection of planets 



on longer time scales. As shown in section 4, the fraction of 
planetary systems from which at least one planet is ejected 
more than doubles in the 10 8 years following a stellar en- 
counter. These later ejections are caused by planet-planet 
scattering, triggered because the fl y-by increased the ec- 
centr icities of the planets (see also [Zakamska fc Tremaind 
l2004h . and/or changed their semi-major axes, leaving the 
system in an unstable configuration. 

(iii) If the intruder in a fly-by is a low-mass star, or even 
a brown dwarf, an exchange reaction may oc cur, similar t o 
what is seen for stellar binaries in clusters (e.g. Heggia ll975l ). 
The intruder may then become bound to the host star, sig- 
nificantly increasing the effect of such fly-bys on planetary 
systems. We have done a large set of scattering simulations 
of encounters between planetary systems and low-mass in- 
truders, and find that such encounters can play an important 
role in the evolution of planetary systems in stellar clusters. 
This shows that a more detailed study of such encounters is 
warranted and we will present that in a forthcoming paper. 

(iv) In general, fly-bys excite the eccentricities of planets, 
both through direct scattering off the intruder, and through 
any subsequent planet-planet scattering. We find that the 
eccentricity distribution of planets which have been scat- 
tered onto orbits tighter than 4.5 au in fly-bys involving 
planetary systems that initially resembles the solar system, 
is similar to that of the observed extrasolar planets (see Fig. 
7 and discussion in section 9.2). It is important to note that 
varying the initial properties of the planetary systems will 
significantly affect the properties of the systems post scatter - 
ing (see also lJuric fc Tremaine!l2008l ; iRavmond et al.ll2009l ). 
We explore this further in section 6 and 7. 

(v) Planet-planet scattering decreases the semi-major 
axes of one or more planets, while leaving others on either 
much wider orbits or ejected. However, planet-planet scat- 
tering in systems resembling the solar system, with gas gi- 
ants initially on orbits of 5 au or wider, can only decrease the 
semi-major axis of the inner planet by roughly a factor of 
two. Hence, the semi-major axes distribution of such post- 
scattering systems does not match the observed distribu- 
tion. Planet-planet scattering models thus need a contribu- 
tion from planetary systems with planets on initially tighter 
orbits than those of the gas giants in the solar system in 
order to reproduce both the observed semi-major axes and 
eccentricity distributions. Hence, disc migration, occurring 
before the stellar encounter, is needed t o explain the ob- 
serve d semi-major axes (see, for example, IMatsumura et all 
l2010l ). 

(vi) During planet-planet scattering, one or more plan- 
ets are successively scattered onto wider and wider orbits, 
until they become unbou nd. The remaining planets are scat- 
tered onto tighter orbits (IScharf fc Menoull2009l ; IVeras et al.l 
2009). During this phase it may be possible to observe plan- 
ets on very wide orbits using imaging techniques. It is im- 
portant to note that this behaviour is independent of what 
mechanism triggered the planet-planet scattering. As such, 
fly-by induced planet-planet scattering will produce planets 
on wide orbits. Furthermore, planets which are captured by 
the intruder stars in fly-bys are often left on rather wide 
orbits, and thus could be observed in imaging surveys. We 
explore both these mechanisms in Section 8. 

In Section 2 we review the encounter rates in 
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young open clusters, measured using iV-body simulations 
(Malmbcrg ct al. 20Q3)- In particular, we discuss the rate of 
fly-bys involving solar-mass stars, as these are most appli- 
cable to solar-system-like planetary systems. 

In Section 3, we describe the properties of planetary sys- 
tems immediately after a fly-by, as measured through both 
numerical simulations and analytical calculations. We mea- 
sure the fraction of systems from which planets are imme- 
diately ejected, and the rate at which planets are captured 
by the intruder star. We also measure how important rela- 
tively distant fly-bys are in terms of moderately exciting the 
inclinations and eccentricities of planetary systems. 

We describe the evolution of the four gas-giants of the 
solar system after a fly-by in Section 4. For example, we 
measure the fraction of such systems in which planet-planet 
scattering is triggered within 10 8 years after the fly-by. We 
also discuss the possibility of predicting the future evolution 
of such systems immediately after the fly-by. 

In Section 5 we measure the rate at which low-mass 
intruders are captured in encounters with planetary systems, 
and discuss the implications such encounters may have on 
the future evolution of the system. 

We proceed to measure the post-fly-by evolution of 
planetary systems slightly different from the solar system. 
In section 7 we go through the effects of fly-bys on systems 
similar to the four gas giants, but where we have changed 
the masses of the planets. In section 8 we simulate the evo- 
lution after a fly-by, of a planetary system consisting of 5 
gas-giants, ge n erate d in simulations of planet formation by 
iLevison et all (|l998l ). 

In section 9 we combine the encounter rates in stellar 
clusters, with the measured effects of fly-bys, to understand 
the effect on the population of planets in clusters of differ- 
ent sizes. We also describe some of the properties of this 
population. 

In section 10, we discuss the implications of encounters 
where the intruding star is in fact an intruding binary. We do 
not consider the effect of s uch encounters here. I n our sim- 
ulations of stellar clusters (jMalmberg et al .1 12007) we had a 
relatively low primordial binary fraction, and as such these 
encounters are not frequent. If we were to increase the pri- 
mordial binary fraction, however, we would see many more 
such encounters, and we discuss how this would affect our 
results. We also, for example, review some of the effects of 
fly-bys on protoplanetary discs, and discuss how changing 
the primordial binary fraction in our clusters can affect our 
results. We summarise our results in section 11. 



2 CLOSE ENCOUNTER RATES IN YOUNG 
STELLAR CLUSTERS 

In this section we review our simulations of stellar clusters, 
and describe the encounter rates of solar-mass stars. 



2.1 Cluster properties 

We have performed a large set of N-body simulations of 
young stellar clusters, with properties representative of the 
open cluster popula tion in the solar neighbourhood (see 
iMalmberg et al.ll2007l . for more details). We used the pub- 
licly available package NBODY6, which is a full force- 



summation direct TV-body code. A complete description of 
the physi cs includ e d and algorithms used in NBODY6 can be 
found in lAarsethl (|2003r ). During the simulations, we mea- 
sured the encounter rates of stars and binaries, focussing 
especially on fly-bys and so-called exchange encounters. In 
our simulations we define a fly-by as when two stars pass 
each other with in a minimum separation of r m i n < 1000 au. 
In an exchange encounter, a single star encounters a binary 
system and replaces one of the stars in the binary. 

In our simulations we consider the evolution of the stel- 
lar clusters after the removal of the primordial gas initially 
present in the cluster. The properties of the clusters were 
chosen so as to represent the observed open clusters in the 
solar neighbourhood. At least 10 per cent of stars form in 
such clusters, while the remaining form either in clusters 
which disperse due to gas removal or in small groups or as- 
sociations (e.g. van den Berghlll98ll; jElmegreen fc Clemens ! 



I l985[ ; iBattinelli fc Capuzzo-Dolcettalll99ll ; lAdams fe Mver; 
120011 ) . Clusters which survive the loss of primordial gas nev- 
ertheless disperse; our simulated clusters have lifetimes of 
between 400 and 1100 million years. 

In this paper, we consider clusters with an initial half- 
mass radius, rh = 0.38 pc and an initial number of stars 
of N = 150, 300, 500, 700 and 1000 respectively. In all these 
simulations, one third of stars were in primordial binaries. 
T he stellar masse s were drawn from the initial mass function 
of iKroupa et~"aH (|l993h with a lower mass of 0.2Mq and 
upper mass of 5Mq. The initial positions of the stars were 
chosen suc h that the i r dist ribution followed the spherically 
symmetric iPlummerl (| 19 1 lT ) model: 



P(r) 



3Af rf 



47rr§ [1 + (r/7-0) 2 ] 5 / 2 : 



(1) 



for total cluster mass M c i. Here, ro is a constant, which 
is connected to th e half-mass radius of the cluster as 
(|Heggie fc Hudl2003ft : r h « 1.305r . 

For each particular cluster type, we performed 10 real- 
isations, performing a new draw from the IMF and gener- 
ating a new set of initial positions and velocities for each 
realisation. 



2.2 Encounter rates 

In Fig. \T\ we show the fractions of initially single stars that 
underwent fly-bys and/or exchange encounters in the cluster 
with N — 700. The results are averaged over 10 realisations. 
The cluster quickly expanded from its initial half-mass ra- 
dius of rh = 0.38 pc to rh ~ 2 pc, whereafter r m i n was 
roughly constant for the remaining lifetime of the cluster. 
From now on, we will refer to this cluster as our reference 
cluster. The circle labelled B contains the stars which spent 
some time in a binary during the cluster's lifetime, the cir- 
cle labelled F contains those stars which underwent a fly-by 
with other stars. Finally, the circle labelled S contains those 
stars which were single at the end of the cluster's lifetime. 
About 6.5 per cent of the initially single stars spent time 
inside a binary system, while about 73.5 per cent underwent 
a at least one fly-by with another star. About 20 per cent 
of the initially single stars did not spend time inside a bi- 
nary or undergo a fly-by. We call such stars singletons. A 
singleton is defined by us as: 
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(i) a star which has not formed in a binary, 

(ii) a star which has not later spent time within a binary 
system, 

(iii) a star which has not suffered close encounters with 
other stars. 

A planetary system orbiting a star that is not a single- 
ton may have had its evolution changed due to interactions 
with other stars. 

In Fig. [2] we show the interaction rates in our reference 
cluster (same cluster as in Fig. [1]), but now only for stars 
with 0.8 Mq < m < 1.2 Mq. As can be seen in the figure, 
the fraction of initially single stars which remain singletons 
throughout the lifetime of the cluster is lower for solar-mass 
stars; only about 15 per cent. The fraction of solar-mass stars 
which spend some time within a binary system is roughly 
similar, while the fraction which have undergone at least one 
fly-by is higher (78 per cent). This increase in the frequency 
of interactions is due in part to the effects of mass segrega- 
tion, which draws the heavier solar-like stars into the denser 
central regions of the cluster. 

In Fig. [3] we plot the distribution of the number of fly- 
bys experienced, ntb, for the initially-single solar-mass stars 
that have undergone fiy-bys in our reference cluster. About 
75 per cent of them have experienced two or more fiy-bys. 
The mean number of fiy-bys per star is in fact four. Stars 
which have undergone several fiy-bys will, on average, have 
closer encounters with other stars than systems only under- 
going a single flyby. The cross section for two stars, having a 
relative velocity at infinity of Voo , to pass within a minimum 
distance r m i n is given by 

o = ^ (l + £) (2) 

where v is the relative velocity of the two stars at closest ap- 
proach in a parabolic encounter i.e. V 2 = 2G(mi+m2)/fmin, 
where mi and ni2 are the masses of the two stars. The sec- 
ond term is due to the attractive gravitational force, and is 
referred to as gravitational focussing. In the regime where 
v <C Voo (as might be the case in galactic nuclei with ex- 
tremely high velocity dispersions), we recover the result, 
a oc r^ lin . However, if v 3> t>oo as will be the case in sys- 
tems with low velocity dispersions, such as young clusters, 
a oc r m i n . Thus a star that has had a number of fiy-bys, no,, 
within 1000 au will on average have passed within a mean 
minimum distance, (r m i n ) min , of: 

V mill 

Lin = lOOOau/nfb. (3) 

As the mean number of fiy-bys in our reference cluster 
is four, we thus find (r m i n ) , = 250 au. 

In Fig. [4]we plot the distribution of when fiy-bys occur, 
£fb (fmin < 1000 au) as a function of cluster age. As can be 
seen in the figure, the evolution is very similar for all clus- 
ters. On average, half of the fiy-bys occur within the first 10 
million years This time-scale is in principle comparable to 
the gas-giant planet formation time-scale, but as our simu- 
lations only treat the post-gas-loss life of clusters, gas-giant 
planets have already formed at the start of our simulations. 



S 




Figure 1. Venn diagram of the stars which were initially single 
in our reference cluster (TV — 700 and T n> initial = 

0.38 pc). The 

numbers are the average result from 10 realisations. The upper 
circle contains stars which were single at the end of the run (S), 
the lower left-hand circle contains the stars which were in a bound 
system (i.e. triple or binary) during the lifetime of the cluster 
(B) and the lower right-hand circle contains the stars which had 
experienced a fly-by (F) (defined as when two stars pass within 
1000 au of each other). 



3 THE PROPERTIES OF PLANETARY 
SYSTEMS IMMEDIATELY AFTER A 
FLY-BY 

We begin here by describing in detail the effects of fiy-bys on 
a planetary system consisting of the four gas-giant planets 
of the solar system (4G). The immediate effect of a fly- by 
depends crucially on the minimum separation between the 
intruder star and the planetary-system-host star during the 
encounter, r m i n . We divide fiy-bys into two different regimes, 
depending on r m i n ; the strong regime (r m i n < 100 au) and 
the the weak regime (100 < r m m < 1000 au). 

To analyse the encounters we use a combination of di- 
rect numerical simulations (in the strong regime) and ana- 
lytic theory (in the weak regime). In all our simulations, we 
set the relative velocity of the two stars, Voo — 1 km/s, which 
is the typical velocity dispersion in young stellar clusters. 

To perform the numerical simulations, we use the 
Bulirsch-Stoer a lgorithm, as im plemented in the MER- 
CURY6 package (|Chamberdll999h . 

3.1 The strong regime: r m i n < 100 au 

Close fiy-bys cause significant changes in the orbits of plan- 
ets. Often, they lead to the immediate ejection of planets 
or an exchange, in which planets become bound to the in- 
truder star. We have performed 6 x 10 5 scattering experi- 
ments of fiy-bys involving the four gas-giant planets of the 
solar system, with r m i n < 100 au and intruder-star mass, 
mi = 0.6, 1.0 and 1.5M . 



© 2010 RAS, MNRAS 000,1111221 



The effects of fly-by s on planetary systems 5 



s 




Figure 2. Venn diagram of the stars which were initially single in 
our reference cluster (TV = 700 and r n j n i t j a i = 0.38 pc) and had 
mass between 0.8Mq and 1.2Mq . The numbers are the average 
result from 10 realisations. The upper circle contains stars which 
were single at the end of the run (S), the lower left-hand circle 
contains the stars which were in a bound system (i.e. triple or 
binary) during the lifetime of the cluster (B) and the lower right- 
hand circle contains the stars which had experienced a fly-by (F) 
(defined as when two stars pass within 1000 au of each other). 




Fly-bys 

Figure 3. The distribution of number of fly-bys, nfl y _b ys , per 
star for the initially single stars with mass 0.8 < m < 1.2Mq 
that experienced fly-bys in our reference cluster (TV = 700 and 
r h, initial = 0-38 pc). The average number of fly-bys per star in 
the cluster is 4. 
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Figure 4. The fraction of all encounters which have occurred as 
function of cluster age, averaged over the 10 realisations for each 
cluster. Only clusters with initial half-mass radii of 0.38 pc are 
included. These clusters have radii of about 2-4 pc during most 
of their lifetime. 




100 



Figure 5. The probability that no planets are ejected or captured 
(thin lines) and that one planet is captured by the intruder star 
(thick lines) immediately after the fly-by as a function of r m ; n 
for encounters involving the four gas giants. The solid lines are 
for encounters where mi = 0.6Mq, the long-dashed lines are for 
encounters where mi = I.OMq and the short-dashed lines are for 
encounters where mj = 1.5Mq. 
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Figure 6. The eccentricities of the the four gas giants of the 
solar system as function of number of fiy-bys, rift,. The minimum 
separation between the host star and the intruder star was chosen 
randomly between 100 and 1000 au. 



In Fig. [5] we plot the probability that four planets re- 
main bound to the host star immediately after the fly-by 
(thin lines) and the probability that one planet is captured 
by the intruder star (thick lines), as a function of r m i n . The 
semi-major axis of the outermost planet, Neptune, is 30 au. 
The probability that all four planets remain bound to the 
intruder star (thin lines) reaches unity at between two and 
three times the semi-major axis of the outermost planet, de- 
pending on the mass of the intruder star. For example, in 
fly-bys with r m i n < 100 au and mi = O.6M0, 15 per cent 
of fly-bys lead to one or more planets becoming unbound 
from the host star. Hence, in fly-bys where r m i n is compara- 
ble in size to the semi-major axes of the planets' orbits, the 
result is often to immediately unbind one or more planets. 
Sometimes planets are ejected from the system and left as 
free-floating planets, but it is not uncommon that planets 
instead are captured by the intruder star, and hence orbit it 
after the fly-by. In the latter case, the most likely outcome 
is that only one planet is captured (thick lines in Fig. [5]). 

3.2 The weak regime: 100 < r min < 1000 au 

In encounters with r m i n > 100 au we can measure the change 
in eccentricities and inclinations of the planets analytically 
(|Heggie fc R asio 1996). Three different cases must be con- 
sidered, depending on the initial eccentricity of the planet 
(circular or eccentric) and on r m i n . 

We have calculated the the average change in the ec- 
centricity of the gas giant planets in the solar system with 
number of fly-bys, nth using equation (|A1|) . (|A2|) and (|A4|) 
from Appendix A. We assume that the intruder is randomly 
oriented with respect to the planets and has mi = I.OMq. 
We plot the result in Fig. [6] A large number of fly-bys is 
required to significantly change the mean eccentricities of 



the planets. Also, as the cross-section for encounters scales 
with r m i n , a star will experience at least one encounter with 
r m in < 100 au for every 10 encounters it experiences with 
100 < r m i n < 1000 au. It is thus unlikely that fly-bys in the 
weak regime trigger planet-planet scatterings in planetary 
system similar to the solar system. They can, however, be 
important as a "heating mechanism" of planetary systems. 
For example, if we assume that the solar system initially 
started out circular and co-planar, we can ask the question 
"How many fly-bys, nth would be required to increase the 
values of the planets' eccentricities and inclinations to their 
current values?" A complication here is that planet-planet 
interactions work to transfer angular momentum between 
planets, continuously changing the eccentricities and incli- 
nations of the planets. Hence, any change in the inclination 
and eccentricities of the planets will propagate throug h the 
system (see, for example. I Zakamska fc Tremainell2004h . 

To estimate the number of fly-bys, no, , required to take 
an initially co-planar and circular version of the solar sys- 
tem and give it similar eccentricities and inclinations to the 
present day solar system, we have chosen to calculate the so- 
called angular momentum deficit of the system, as a function 
of n fb . 

The angular momentum deficit (AMD) in a system is 
the additional angular momentum present in it due to the 
non-circularity and non-co-planarity of the planetary orbits. 
It can be thought of as the am ount of non- linearity present 
in the system. It is defined as (|Laskarll2OQ0l ): 

AMD = ^ A k (l — -y/l — e£ cos i k ) (4) 

k=l 

where ek and ik are the eccentricity and inclination of the 
fc th planet, and A = m^^/ua-^, where and are the 
mass and semi-major axis of the fcth planet. Starting out 
with a co-planar system with planets on initially circular 
orbits we find that 15 fly-bys with solar-mass intruders with 
100 < r m i n < 1000, is enough to give the system its present 
AMD value. 



4 POST-ENCOUNTER EVOLUTION OF 
PLANETARY SYSTEMS 

In this section we consider the evolution of the planetary sys- 
tems on the longer time-scale after fly-bys. In multi-planet 
systems, planet-planet interactions will cause evolution of 
the planets' orbital elements. The effect of most fly-bys (ex- 
cept the very tight ones) is to only change the angular mo- 
menta of the planets, with the largest change for the outer- 
most planet. While this change may not immediately trigger 
planet-planet scattering, redistribution of angular momen- 
tum between the planets (e.g. planet-planet interactions) 
may cause the eccentricities of one or more planets to grow 
over time, triggering planet-planet scattering. 

4.1 The long-term evolution of the four gas giants 
of the solar system after a fly-by 

We begin here by describing the long-term evolution of plan- 
etary systems consisting of the four gas giants of the solar 
system that have undergone fly-bys with r m i n < 100 au. 
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Table 1 . Overview of the evolution of the four gas giants of the solar system after a fly- by, divided into different categories depending 
on the state of the system immediately after the fly-by. In column 2 we list the categories and in column 3 we list the fraction of all 
systems which belongs to a given category. In column 4, 5, 6 and 7 we divide the systems from each category into four bins, depending 
on at what time after the fly-by the first planet was ejected from the system, t c j ect . In column 4 are systems with t e j ec t < 10 Myr, in 
column 5 systems with teject < 30 Myr, in column 6 systems with t e ject < 100 Myr and in column system are systems from which no 
planet was ejected in 10 8 years. Its worth pointing out, that for systems in which a planet is ejected immediately in the fly-by, column 4 
and 5 are by definition 0, as all such planets have t e ject = 0. Furthermore, the fraction of systems which falls into the categories "Kozai 
immediately after fly-by", "AMD unstable immediately after fly-by" and Kozai before close encounter" are in general rather small the 
uncertainty in these values is rather large, especially in column 3-6. For example, in the simulations with rrtj = 0.6Mq, only 1.8 per cent 
of systems, or 18 simulations, fall into the category AMD unstable immediately after fly-by", with 3 systems that have i c jcct < 10 7 years 
(column 3), 2 systems with 10 7 < t c j C ct < 3 X 10 7 (column 4), 6 systems with 3 X 10 7 < teject < 10 s (column 4) and finally 7 systems 
which are AMD unstable immediately after the fly- by, but from which no planets are ejected within 100 Myr after the fly- by (column 6). 



mi Situation immediately t c jcct (Myr) No planet ejection 

after flyby 10 30 100 after 100 Myr 



Immediate ejection: 


0. 


.150 


1. 


.00 


0. 


00 


0. 


.00 


0. 


00 


livrut" r"T"(~it;t;i n cr iminoniatolv cifrrM" fliz— rv\/" 





.262 





.62 


0. 


18 





.20 


0. 


00 


Kozai immediately after fly-by: 





.010 





.10 


0. 


00 


0. 


.10 


0. 


.50 


AMD unstable immediately after fly-by: 





.018 





.17 


0. 


11 


0. 


.33 


0. 


39 


Kozai before close encounter: 


0. 


.004 





.25 


0. 


00 


0. 


.00 


0. 


75 


Other (unstable): 





.037 





.10 


0. 


22 


0. 


.68 


0. 


00 


Othcr(stablc): 





.519 





.00 


0. 


00 





.00 


1. 


00 


Immediate ejection: 





.251 


1 


.00 


0. 


00 


0. 


.00 


0. 


00 


Orbit crossing immediately after fly-by: 





.289 





.58 


0. 


20 


0. 


.22 


0. 


00 


Kozai immediately after fly-by: 





.012 


0. 


.25 


0. 


08 


0. 


.17 


0. 


50 


AMD unstable immediately after fly-by: 





.015 





.00 


0. 


27 





.27 


0. 


46 


Kozai before close encounter: 





.001 


1. 


.00 


0. 


00 


0. 


.00 


0. 


00 


Other (unstable): 





.039 





.00 


0. 


15 


0. 


.85 


0. 


00 


Other(stablc): 





.393 





.00 


0. 


00 





.00 


1. 


00 


Immediate ejection: 





.308 


1 


.00 


0. 


00 


0. 


.00 


0. 


00 


Orbit crossing immediately after fly-by: 





.304 





.58 


0. 


20 


0. 


.22 


0. 


00 


Kozai immediately after fly-by: 





.022 





.23 


0. 


14 


0. 


.08 


0. 


.55 


AMD unstable immediately after fly-by: 





.031 





.06 


0. 


19 





.35 


0. 


40 


Kozai before close encounter: 


0. 


.005 





.10 


0. 


00 


0. 


.00 


0. 


60 


Other (unstable): 





.050 





.16 


0. 


30 


0. 


.54 


0. 


00 


Othcr(stablc): 





.280 





.00 


0. 


00 





.00 


1. 


00 



We have performed three sets of simulations, with mi = 
0.6, 1.0 and 1.5M Q respectively). For each set we perform 
1000 realisations, between which we randomly vary the ori- 
entation of the intruder star's orbit. We vary r m i n linearly 
between and 100 au, which is consistent with the encounter 
rates we measure in our simulations of young stellar clus- 
ters (Section [2]). In each simulation we measure the time 
until the first close encounter between two planets occur, 
tencountcr, and the time until one planet is ejected teject- A 
close encounter is defined as when two planets pass within 
one mutual Hill radii, Ru, of each other, where: 



mi = I.5M0). The fraction of systems in which a planet 
is ejected immediately after the fly- by increase with mi. 
For mi = 0.6, 1.0 and 1.5Mq 15, 25 and 31 per cent of 
systems suffer an immediate ejection. The fraction of sys- 
tems from which at least one planet has been ejected then 
grows with time, reaching values of 0.47 (mi = O.6M0), 0.59 
(mi = l.OA/ ) and 0.69 (mi = 1.5M©) at 10 s years after 
the fiy-by. Hence, planet-planet scattering, triggered due to 
the fly-by, more than doubles the fraction of systems from 
which at least one planet is ejected in 10 s years for the 4G 
system. 



/ \ 1/3 

Hi + a 2 / mi+m 2 \ 

where mi and m 2 are the the masses of the planets, ai and 
02 are their semi-major axes and Mh is the mass of the host 
star. 

In Fig. [T^i we plot the fraction of systems from which 
at least one planet has been ejected, /(teject), as a func- 
tion of time after the fiy-by. The three different lines corre- 
spond to different intruder star masses (solid line: mi = 
O.6M0, dashed line: mi = I.OJW0 and dash-dotted line: 



In Fig. 0) we plot the fraction of systems in which two 
planets have undergone a close encounter with each other, 
/^encounter). We define a close encounter as when two plan- 
ets pass within one mutual Hill radius of each other (see 
equation[5]). We find that /(tencounter) is always significantly 
larger than /(t e ject), implying that the ejection of a planet, 
unless it occurs immediately after the fly-by, is preceded by 
a phase of planet-planet scatterings. As /(tencounter) is larger 
than /(teject) at t = 10 8 years we can also predict that the 
fraction of systems from which at least one planet has been 
ejected will continue to increase for t > 10 s years. 
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Figure 7. The evolution of the four gas giants of the solar system, 
after a fly-by with r m i n < 100 au. The periastron passage of the 
intruder star occurred at t = 0. In panel a) we plot the fraction 
of systems, /(t e ject)i from which at least one planet has been 
ejected within, f j ec t years after the fly-by and in panel b) we 
plot the fraction of systems, /(tcncountcr), in which at least two 
planets have passed within one mutual Hill radius of each other 
as a function of time, t cnC ountcr after the fly-by. 



4.2 Predicting the long-term effects of fiy-bys 

While one in principle could simulate the long-term evolu- 
tion of planetary systems which have undergone fiy-bys for 
a large range of planetary system architectures and fly-by 
properties, it would be greatly beneficial if one could predict 
the future evolution immediately after the fly-by. 



We have broken down the post-encounter systems from 
our simulations of the four gas giants into seven different 
categories, depending both on their properties immediately 
after the fly-by and on their long-time evolution. The sys- 
tems are put into only one category; the first one in the 
following lists whose criteria they fulfil. 

(i) Immediate ejection: planetary systems in which at 
least one planet became unbound from the host star im- 
mediately after the fly-by, 

(ii) Orbit crossing immediately after By-by: planetary 
systems in which the orbits of at least two planets cross 
immediately after the fly-by, 

(iii) Kozai immediately after By-by: planetary systems in 
which the mutual inclination of at least two planets is be- 
tween 39.2° and 140.8° immediately after the fly-by, 

(iv) AMD unstable immediately after By-by: planetary 
systems in which the angular momentum deficit (AMD) im- 
mediately after the encounter is such that the orbits of two 
planets may cross, 

(v) Kozai before close encounter: planetary systems that 
become unstable in 10 s years in which the mutual inclination 
between at least two planets attains a value between 39.2° 
and 140.8° before the onset of planet-planet scattering, 

(vi) Other (unstable): planetary systems which do not ful- 
fil any other criteria, but from which at least one planet is 
ejected within 10 8 years after the fly-by, 

(vii) Other (stable): planetary systems which do not fulfil 
any of the above criteria and are stable for at least 10 s years 
after the fly-by. 

These criteria are described in further detail below. 

If the mutual inclination between two planets is b etween 
39.2° and 140.8° the so-called Kozai mechanism l|Kozail 
I1962T I operates. The effect of the Kozai mechanism is to 
cause the eccentricity and inclination of planets to oscillate, 
triggering planet-planet scatterings in multi-plan et systems 
l|Malmberg et al.ll2007l ; iMalmberg fc Davies|[2009l ). 

While the intruder star in a fly-by will not by itself 
cause the Kozai mechanism to operate, it may leave the 
system in a state where the mutual inclination between two 
or more planets is between 39.2° and 140.8°. The Kozai 
mechanism will th en operate, potentially triggering planet- 
planet scatterings (|Nagasawa et al]|200ct ). 

It is also possible that planet-planet interactions, which 
change both the eccentricities and inclinations of planets 
with time, increase the mutual inclination between two plan- 
ets such that it reaches values between 39.2° and 140.8° even 
if this is not the case immediately after the fly-by. For sys- 
tems which do not fulfil the criteria of category (i)-(iv), but 
in which planet-planet scatterings occur within 10 8 years af- 
ter the fly-by, we check whether this happens, and if so, the 
systems are put into category (v). 

The angular momentum deficit (AMD, see Equation [4} 
in a system is the additional angular momentum present 
in it due to the non-circularity and non-co-planarity of the 
planetary orbits. It can be thought of as the amount of non- 
linearity present in the system. We use the AMD to cal- 
culate whether the planetary system may become unstable, 
by checking if any pair of planets can achieve high enough 
eccentricities for their respective orbits to cross, if all the 
availa ble AMD is put into these two planets (see lLaskarl 
2000, for more details). 
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In column 3 of Table [T] we list the fraction of systems in 
our simulations of fly-bys belonging to each of the categories 
listed above. In columns 4-6 we list the fraction of systems 
from which at least one planet is ejected within 10 7 , 3 x 10 7 
and 10 8 years respectively after the fly-by. Finally, in column 
7, we list the fraction of systems from which no planets are 
ejected in 10 s years after the encounter. 

We can now in more detail analyse the post-encounter 
behaviour of the planetary systems. For example, in encoun- 
ters with an intruder of mass O.GMq, at least one planet is 
ejected immediately after the fly-by in 15 per cent of the 
encounters. An additional 26.2 per cent of all systems are 
left in a state where the orbits of at least two planets cross. 
Thus, in total, 41.2 per cent of systems fall into categories 
(i) and (ii). A total of 47 per cent of systems eject at least 
one planet in 10 s years after the fly-by. Hence, about 88 per 
cent of the systems which become unstable in 10 8 years after 
the fly- by, fall into categories (i) and (ii). This behaviour is 
similar also for fly-bys with intruders of mass I.OMq and 
1.5M . 

It may seem surprising that some systems left with 
orbit-crossing planets do not eject any planets until 100 Myr 
after the flyby encounters. In many cases, planets left on ec- 
centric orbits are also rather inclined to the original plane 
of the unperturbed planetary system, so even if orbits cross, 
close encounters between two planets may be extremely rare 
events. Only a small fraction of systems that later become 
unstable are left in such a state that the Kozai mechanism 
would operate or that they are AMD unstable. About 4-5 
per cent of all systems fall into the category unstable (other). 

Given that the vast majority of systems which eject 
planets fall into categories (i) and (ii), we may use the list 
of post-flyby categories to predict which systems will eject 
planets. It is very useful that we are able to predict the post- 
flyby evolution of a system containing the four gas-giants, 
as it reduces the computational effort needed to understand 
the effect of fly-bys. As we shall in section 6 however, pre- 
dicting the post-fly-by evolution of planetary systems which 
are inherently less stable than the four gas giants is a more 
difficult task. 



5 CAPTURE OF LOW-MASS INTRUDERS 

In this section we consider the effects of encounters involv- 
ing brown dwarfs and free-floating planets. Such encoun- 
ters are different from fly-bys by stellar-mass intruders, as 
the intruder may be captured. These sorts of encounters 
are rather similar in nature to encounters between single 
stars and hard stellar binaries - that is, encounters where 
the binding energy of the binary exceeds the kinetic energy 
of the incoming third star. Ultimately, such encounters re- 
sult in the ejection of one of the three stars, leaving the 
binary more tightly bound. In the case of encounters with 
a star possessing a planetary system, the encounter can be 
somewhat complicated given the larger number of bodies in- 
volved. If the incoming planet is captured, even temporarily, 
close encounters between it and other planets in the system 
can occur, leading to the ejection of planets, or at least leav- 
ing them on altered orbits. 

Roughly speaking, a transient 3-body system can only 
form if the binding energy of the planetary system is larger 



than the kinetic energy of the intruder star. We can write 
the critical velocity, i> cr it, for which the total energy of the 
intruder and planetary systems equals zero, as: 



/ Gm H m p 
V Map 



1/2 



(6) 



where mn is the mass of the planet host star, m p is the 
mass of the planet, mi is the mass of the intruder, [i — 
(toh + m p )mi/(mH + m p + mi) is the reduced mass of the 
planet-intruder system and a p is the semi-ma jor axis of the 
plane t. If Voo < v cr u a capture can occur (|Fregeau et al.l 
2006). For example, from equation [5] we can calculate that if 
Voo — 1-0 km/s, m p = Mj and ap = 30 au, then mi < 30mj. 

We have calculated numerically the capture cross- 
section for sub-stellar mass intruders (mi = 3, 10 and 30m,j) 
in an encounter with a planetary system consisting of a 
Jupiter-mass planet orbiting a solar-mass star at 30 au us - 
ing the TV-body code Starlab (|Portegies Zwart et al.ll200ll ). 
In Table [2] we list the measured cross-sections, cr, scaled as 
the dimensionless quantity: 



^scaled 



(7) 



where v is defined as v = ^oo/ucrit- 

We divide the measured cross-sections into several cat- 
egories, depending on the nature of the encounter. If a tran- 
sient 3-body system forms, the encounter is categorised as 
a resonant interaction, which results in either the planet or 
the intruder becoming unbound. It is also possible that the 
outcome of the encounter is a prompt exchange, where the 
intruder can take the place of either the planet or the host 
star. Finally, if v > 1 the total energy of the planet-intruder 
system is greater than zero, and an ionisation may occur, 
leaving all three objects unbound. 

The measured cross-sections in Table [2] are small if 
Voo ~ 1 km/s. We have also calculated the cross-sections 
for Voo = 0.5 km/s for which the cross-sections are an order 
of magnitude larger. 

We have used the cross-sections in Table [2] to calcu- 
late the fraction of solar-mass stars hosting a Jupiter-mass 
planet with semi-major axis equal to 30 au, that are involved 
in resonant encounters or direct exchanges with sub-stellar 
mass objects in our reference cluster (iV = 700 stars and 
rh.initial = 0.38 pc). 

We use the log-normal IMF o f lChabrierj (|200ll . |2002| ) to 
extend the stellar population in our reference cluster down 
to one Jupiter-mass objects. From our iV-body simulations, 
we find that for encounters involving solar-mass stars and 
intruder stars with mi < IMq, the number of encounters 
that a solar-mass star experiences with an intruder of mass 
mi, is proportional to the number of stars with mass mi 
in the clusters. We assume that this relation holds also for 
sub-stellar mass objects and that the encounter rate for sub- 
stellar mass objects scales linearly with r m i n , as is the case 
for stellar-mass intruders. 

We use three different mass ranges of the intruders in 
order to account for how the cross-sections change with mass 
of the intruder: lMj < m < 5Mj for 3Mj intruders, 5Mj < 
m < 15Afj for lOMj intruders and 15Mj < m < 45Mj 
for 30Mj intruders. We list the fractions of solar-mass stars 
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Table 2. The calculated cross-sections (column 2-4) that encounters between a sub-stellar mass object and a Jupiter-mass planet orbiting 
a solar-mass star with semi-major axis of 30 au, leads to the formation of a transient 3-body system, called a resonant interaction, or 
to a prompt exchange. For intruder mass mj = 30Mj and Uoo = 1.0 km/s, the encounter sometimes lead to ionisation, where all the 
three objects become unbound. We have calculated the cross-section for three different values of m\\ 3, 10 and 30Mj, and two different 
values of Doc,: 0.5 km/s (top panel) and 1.0 km/s (lower panel). In column 5 we list the fraction of such encounters, / c which occur 
in our reference cluster (jV = 700 stars and ; = 1000 au). The cross-sections, a has been scaled to give the dimensionless quantity 
("■scaled = l,2 °'/ 7ra p)i where a p is the semi- major axis of the planet and v = «oo/"crit- 



Type mi = 3mj mi = 10m.j mi = 30m j f e 

Voo = 0.5 km/s 



Resonant interaction, intruder ejected 


8.0±1.3 XlO 


-2 


2.1±0.3 xlO" 


-l 


7.3 ±0.5 XlO" 


-i 


4.1 x 10" 


-2 


Resonant interaction, planet ejected 


1.0±0.4 xlO" 


-2 


0.7±0.2 xl0~ 


-l 


5.6 ±0.5 X10" 


-l 


3.2 x 10- 


-2 


Prompt exchange, planet ejected 


0.0 




0.0 




1.5 ±0.3 X10" 


-l 


1.3 x 10" 


-2 


Prompt exchange, host star ejected 


0.0 




0.0 




0.0 









Ionisation, system unbound 


0.0 




0.0 




0.0 









V 


0.159 




0.291 




0.504 








Uoo = 1.0 km/s 


















Resonant interaction, intruder ejected 


9.8±2.0 xlO 


-3 


1.8±0.6 xlO" 


-2 


0.0 




2.0 x 10- 


-3 


Resonant interaction, planet ejected 


3.0±1.0 xlO" 


-3 


1.4±0.5 xicr 


2 


0.0 




1.5 x 10~ 


-3 


Prompt exchange, planet ejected 


0.0 




8.8±5.0 XlO" 


-4 


1.4±0.2 Xl0" 


-l 


6.5 x 10" 


-3 


Prompt exchange, host star ejected 


0.0 




0.0 




1.5±0.8 X10" 


-2 


2.0 x 10- 


-3 


Ionisation, system unbound 


0.0 




0.0 




5.2±5.0 X10" 


3 


1.2 x 10- 


-3 


V 


0.318 




0.582 




1.007 









which have been involved in these encounters, f c , in column 
5 of Table [2] The exact lower-mass cut-off in the IMF is 
not yet well constrained from observations, a nd may vary 
betw een different clusters (see, for example, IScholz et al.l 
120091 . and references therein) . However, we find that the frac- 
tion of stars involved in these encounters does not change 
significantly if we, for example, set the cut-off at 5Mj in- 
stead. The cause of this is two-fold. Firstly, the capture- 
cross-sections increase rather strongly with mi, and hence 
are smaller for intruder of mass 3Mj than for the higher 
mass intruders. Secondly, given the IMF used, the number 
of objects in the range 1M.j < m < 5Mj is significantly 
smaller than the number of objects with, for example, mass 
15Mj < m < 45Mj. 



It is important to note that we have here assumed that 
all solar-mass stars host planetary systems. Depending on 
the velocity dispersion in the cluster, the fraction of planet- 
hosting stars involved in resonant or prompt exchange en- 
counters with sub-stellar mass objects ranges from a few per 
cent (lice = 0.5 km/s) down to 1CP 3 (uoo = 1.0 km/s). These 
fractions are less than the fraction of four gas-giant systems 
from which planets are ejected due to fly-bys with stellar- 
mass intruders, which is about 11 per cent in our reference 
cluster (see Section 9). The capture of a sub-stellar mass 
object is likely to significantly change the orbits of planets 
when it occurs. Hence, depending on the velocity dispersion 
in clusters, several per cent of planetary systems in stel- 
lar clusters could be significantly changed by such events. 
We will be carrying out more detailed simulations of rates 
at which sub-stellar-mass objects are involved in encounters 
with planetary systems, and how they affect the evolution 
of the system, in a forthcoming paper. 



6 CHANGING THE MASS DISTRIBUTION OF 
THE PLANETS 

In this section we consider the effects of fly-by encounters 
involving planetary systems containing planets with masses 
different from the four gas giants, but with the same initial 
orbits. We have used two different mass-distributions: the 
four Jupiter system (4J), where we set the mass of all the 
gas giants equal to 1 Mj and the geometric mean (GM), 
where we set the masses of the gas giants, m* as: 



m* = VMjrm. (8) 

where m\ is the mass of the i th planet in the normal solar 
system. We call systems in which the masses of the gas giants 
are all rather similar democratic. 

We have simulated the effects of fly-bys on these sys- 
tems numerically, with r m i n < 100 au and mi = 0.6Mq. 
We performed 1000 realisations for each system, varying the 
initial orientation of the intruder star, as well as the initial 
positions of the planets in their respective orbits. In Fig. 
[8] we plot the fraction of systems from which at least one 
planet has been ejected as a function of time after the fly-by 
for the 4G (solid line), the 4J (dotted line), and the GM 
(long-dashed line). Both the 4G and GM are stable if left 
alone. However, in the 4J-system, planet-planet scatterings 
are sometimes triggered without external perturbation. We 
include, in Fig. the fraction of systems from which at least 
one planet has been ejected as a function of time in the 4J 
without a fly- by (short-dashed line). 

In the case of the 4J, about 8 per cent of systems be- 
come unstable, and eject at least one planet in 10 8 years 
without a fly-by. However, with a fly-by, the rate at which 
the 4J-systems become unstable and eject planets is greatly 
enhanced. About 73 per cent of the 4J-systems which have 
suffered a fly-by with r m i n < 100 au and mi < O.6A/0 be- 
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Table 3. Overview of the evolution of the four gas giants after a fly-by, where now either all four planets have the mass of Jupiter (4J) 
or the planetary masses are each the geometric mean of Jupiter and their original mass (GM). We have divided the systems into different 
categories, depending on the state of the system immediately after the fly-by. In column 1 we list the categories and in column 2 we list 
the fraction of all systems which belongs to a given category. In column 3, 4, 5 and 6 we divide the systems from each category into 
four bins, depending on at what time after the fly-by the first planet was ejected from the system, t G j cc t - In column 3 are systems with 
tcjcct < 10 Myr, in column 4 systems with t e j Gct < 30 Myr, in column 5 systems with t e j ect < 100 Myr and in column 6 are systems 
from which no planet was ejected in 10 s years. Its worth pointing out, that for systems in which a planet is ejected immediately in the 
fly-by, column 4 and 5 are by definition 0, as all such planets have t e j ect = 0. Furthermore, the fraction of systems which falls into the 
categories "Kozai immediately after fly-by" , "AMD unstable immediately after fly-by" and Kozai before close encounter" are in general 
rather small the uncertainty in these values is rather large, especially in column 3-6. For example, in the case of the geometric mean 
system, only 1.6 per cent of systems, or 16 simulations, fall into the category AMD unstable immediately after fly-by", with 3 systems 
that have t e ject 

< 10 7 years (column 3), systems with 10 7 < t e j ect < 3 X 10 7 (column 4), 6 systems with 3 X 10 7 < t 

eject *^ 10^ (column 

4) and finally 7 systems which are AMD unstable immediately after the fly-by, but from which no planets are ejected within 100 Myr 
after the fly- by (column 6). 
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Immediate Ejection: 


0. 


153 


1. 


00 





.00 


0. 


00 


0. 


.00 


Orbit crossing immediately after fly-by: 


0. 


263 


1. 


00 





.00 


0. 


00 


0. 


.00 


Kozai immediately after fly-by: 


o. 


017 


0. 


80 


0. 


.13 


0. 


00 





.07 


AMD unstable immediately after fly-by: 


o. 


045 


0. 


80 


0. 


.03 


0. 


08 





.10 


Kozai before close encounter: 


0. 


000 


















Other (unstable): 


0. 


264 


0. 


60 





.16 


0. 


24 


0. 


.00 


Other (stable): 


0. 


258 


0. 


00 





.00 


0. 


00 


1 


00 



etary system with different number of planets and different 
orbital elements than the gas giants in the solar system. 

We have used one of the syste ms generated insim- 
ulations of planetary formation by iLevison et alj (|l998h . 
namely the system generated in run 10 of series A (hence- 
forth called RlOa). The orbital properties of the planets in 
this system can be found in table |4j and was chosen as a 
system to be broadly similar to the solar system, with gas 
giants on orbits between about 5 and 30 au. This system 
is stable for at least 10 9 years. It contains five planets, and 
is slightly more tightly packed than the solar system. The 
masses of the planets are slightly more democratic than that 
of the four gas giants of the solar system. We have per- 
formed three sets of fly-by simulations, with r m i n < 100 au 
and mi = 0.6Mq,1.0Mq and 1.5Mq respectively. For each 
set, we have done 200 realisations, varying the initial orien- 
tation and r m i n of the intruder star. We plot the fraction of 
systems from which at least one planet has been ejected in 
Fig. [9] As can be seen from this figure, this system is vul- 
nerable to fly-bys. For example, fly-bys with mi = 0.6Mq 
cause 64 per cent of these systems to undergo planet-planet 
scattering and eject at least one planet in 10 s years, which 
is a number somewhat larger than seen previously for the 
four gas giants in the solar system. This is not surprising, 
as the system is a more packed and more democratic. 

Unsurprisingly, we find that it is the two outer (and 
lowest mass) planets that are ejected most often. However, 



come unstable and eject at least one planet within 10 s years 
of the fly-by. For the GM systems, about 56 per cent of 
systems eject at least one planet in 10 8 years after a fly-by. 

We have carried out the same analysis of the GM and 
the 4J-systems as we did for the 4G in section 14.21 (see ta- 
ble [3]) As can be seen, the fraction of systems from which 
a planet is ejected immediately after the fly-by, and the 
fraction of systems in which the orbits of planets are found 
to cross, is similar to that which we find from our simula- 
tions of the four gas giants. However, the fraction of systems 
that become unstable, but do not fit into any specific cate- 
gory, is much larger. About 12 per cent of GM-systems and 
about 26 per cent of the 4J-systems fall into the category 
Other (unstable). Hence, we could not predict that these sys- 
tems would become unstable by analysing the orbits of the 
planets immediately after the fly-by. 



7 VARYING THE ARCHITECTURE OF THE 
PLANETARY SYSTEMS 

So far, we have only considered the effects of fly-bys on plan- 
etary systems derived from our own solar system. However, 
as the population of planets around stars is likely to also 
include systems with a different architecture to our solar 
system, we have also simulated the effect of fly-bys on a plan- 
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Table 4. The orbital elements of the planetary systems generated in run 10, series A of iLevison et al ] jl99Sh (RlOa). The columns 
are mass m (in Jupiter Masses, Mj), semi-major axis a (AU), eccentricity e, inclination i, argument of perihelion cj, longitude of the 
ascending node Q, and mean anomaly M. 





m (Afj) 


a (AU) 


e 


i 


(J 


n 


M 


Planet 1 


0.880 


4.92321 


0.02742 


1.49983 


12.14303 


108.14294 


90.04206 


Planet 2 


0.618 


8.34403 


0.00684 


1.18543 


40.16036 


279.69692 


28.16437 


Planet 3 


0.140 


11.78578 


0.03653 


2.79969 


174.64686 


0.49883 


132.94036 


Planet 4 


0.060 


18.69599 


0.03278 


3.28732 


342.47166 


244.48202 


133.67521 


Planet 5 


0.041 


32.59675 


0.02001 


0.98343 


25.15793 


287.33401 


309.54848 




Figure 8. The evolution of the four gas giants of the solar system 
(4G) after a fly-by, the four Jupiters (4J), both if left on its own 
and after a fly-by and the evolution of a systems where the mass of 
the gas giants in the solar systems was set as the geometric mean 
(GM) between the 4G and the 4J after a fly-by. In all simulations 
we varied r m i n linearly between and 100 au and the mass of 
the intruder star was mi = 0.6Mq. The periastron passage of the 
intruder star occurred at t = 0. We plot the fraction of systems, 
/(teject)i from which at least one planet has been ejected within, 
ieject years after the fly-by 



at least three planets are ejected in about 35 per cent of 
fly-bys with mi = O.GMq. As seen with earlier simulations 
of other planetary systems, the planets left behind have a 
wide range of eccentricities. The RlOa system used here, is 
thus another example of a stable planetary system which 
is often destabilised via encounters in stellar clusters. The 
solar system is not unique; there could be a large family of 
planetary architectures that follow this behaviour. 



8 PLANETS ON VERY WIDE ORBITS 

A number of planets have been observed on very wide orbits 
around ma in-sequence dwarf stars by direct im aging tech- 
niques (e.g. lKalas et al.ll2008l ; lMarois et al.ll200Sl ). These im- 



Figure 9. The evolution of the plan etary system produced in run 
10 of series A in lLevison et all lll998T) after a fly-by. The properties 
of the intruder star was r m ; n < 100 au and mi = O.GMq (solid 
lines), I.OMq (dashed line) and 1.5Mq (dashed-dotted line). The 
periastron passage of the intruder star occurred at t = 0. 



aged planets have projected separations to their respective 
host stars ranging from 25 to about 300 au. The host stars 
of these planets are A-stars with masses of 1.5 — 2.0AfQ. 

It is hard to form planets at such large distances from 
the host star in the core-accretion scenario, as the surface 
density of the disc is too low. In this section we explore how 
fly-bys can put planets on very wide orbits. 

8.1 Planet-planet scatterings 

Gas-giant planets which form at smaller orbital radii, can 
be put on very wide o rbits by planet-p lanet scattering 
|Scharf fc Menoul 120091 ; IVeras et all l2009h . Planet-planet 
scattering often lead to the ejection of one or more plan- 
ets. However, in general planets are not ejected in a sin- 
gle scattering event, but rather suffer a series of encounters 
which put them on successively wider orbits, until they fi- 
nally become unbound. This process can take several 10 7 
years, during which the planets can be observed on very 
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semi— major axis (au) 



Figure 10. The semi-major axes, a and eccentricities, e of all 
planets which are bound to the host star 10 8 years after the fly- 
by in our simulations of fly-bys involving the four gas-giants of 
the solar system with r m ; n < 100 au and mj = 0.6Mq . 



also likely that at least t wo of the planets are locked in a 2:1 
mean motion resonance (|Fabrvckv fc Murrav-Ciavll201Ch . 

In Fig. [10] we show a snapshot of the semi-major axes 
and eccentricities of the four gas-giant planets in the solar 
system, 10 s years after a fly-by with r m i n < 100 au and 
mi = 0.6A/q. At this time, about five per cent of systems 
contain at least one planet with semi-major axis larger than 
100 au, which are planets that could be observed in imaging 
surveys. 

In Fig. [TTJ we plot the fraction of systems that have at 
least one planet with semi-major axis larger than 100 au 
as a function of time after a fly-by involving the 4G, with 
f'min < 100 au and mi = Q.QMq. As can be seen, the fraction 
of planets on wide orbits varies with time. It reach its peak 
value (between 0.07 and 0.11 depending on the mass of the 
intruder star) in only a few 10 years after the fly-by and 
then slowly decrease with time. 

From Fig. [TTJ one can see that already at t = 10 4 years 
there are planets on wide orbits. These planets were not 
put on such wide orbits by planet-planet scattering. Instead, 
their wider orbits were caused by the planets strongly inter- 
acting with the intruder star in close fly-bys. The fraction 
of systems with planets on wider orbits is largest for fly- 
bys with mi = O.6M0. The cause is that in these fly-bys, 
the perturbation on the planetary orbits is weaker than for 
higher mass intruders. This decreases the fraction of sys- 
tems from which planets are ejected promptly. However, it 
increases the fraction of planets which are almost ejected, 
and hence are put on wide orbits. 



— m,=0.6M s 

— - m,= 1.0M E 

— m,= 1.5M s 

0.15 - 




10" 10 s 10" 10' 10 s 

t 



Figure 11. The fraction of systems, /p, containing planets with 
semi-major axes greater than 100 au, as a function of time, t, 
after the fly-by, in our simulations of fly-bys involving the four 
gas giants of the solar system with r m ; n < 100 au. 



8.2 Captured planets 

In tight fly-bys, where the minimum separation between the 
host star and the intruder is only a factor two or three larger 
than the semi-major axis of the outermost planet, one or 
more planets may become bound to the intruder star dur- 
ing the encounter. In Fig[S]we plot the probability that one 
planet is captured by the intruder star in encounters involv- 
ing the four gas giants of the solar system, where r m i n < 100 
au (thick lines). As an example, in our reference cluster, the 
fraction of stars in the mass-range 1.3 < m < 1.7Mq that 
would capture one planet in fly-bys is 0.08 x / p , where / p is 
the fraction of stars in the cluster hosting planetary systems 
similar to the 4G. 

Planets which are captured during a fly-by are in gen- 
eral left on rather wide and moderately eccentric orbits. In 
Fig. \T%\ we plot the semi-major axes and eccentricities of 
the planets captured by a 1.5Mq intruder. About 40 per 
cent have semi-major axes larger than 100 au. If a fly-by oc- 
curs when the planets are newly formed, and hence are still 
contracting, the luminosity of the planets, combined with 
their large separations to the host star, should make them 
detectable in imaging surveys. 



wide orbits. It is important to note that while planet-planet 
scatterings can explain the orbits of single planets on wide 
and moderately eccentric orbits (like e.g. Formalhaut b), it 
is very unlikely that systems like HR8799 formed this way. 
HR8799 consists of t hree planets which all appear to be on 
rather circular orbits (jMarois et al.ll2008r ). Furthermore, it is 



9 PROPERTIES OF PLANETARY SYSTEMS 
PRODUCED IN STELLAR CLUSTERS 

In this section we describe how the rate at which planets are 
ejected from solar-system like planetary systems vary with 
the properties of the stellar cluster. Of course, what we are 
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Figure 12. The semi-major axes, a, and eccentricities, e, of 
planets that were captured by the intruder star in fly-bys with 
r min < 100 au, involving the four gas-giant planets of the Solar 
system and a intruder star with mass 1.5Mq. We here only in- 
clude planets in single-planetary systems, i.e. in which only one 
planet was captured by the intruder star. Most of the planets are 
on rather wide orbits, with about 40 per cent having a > 100 au. 
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Figure 13. The fraction of planet-hosting systems from which 
at least one planet is ejected in 10 8 years due to fly-bys. The 
fractions were calculated assuming a single value of the intruder 
star mass of rrtj = 0.6Mq, with r m ; n < 100 au. The error bars 
are the standard deviation from the mean of / c j cc t as measured 
for the 4G. In the figure we have, for clarity, offset N slightly to 
the right for all systems except the 4G. 



Table 5. The fraction of planetary systems of the four gas gi- 
ants (4G), four Jupiters (4J), the geometric mean of the two 
(GM) and the plane tary system produced in run 10 of series A in 
iLevison et al.l (l998) (RlOa) from which planets are immediately 
ejected, /ion, captured by the intruder star, / ca p, and the total 
fraction of systems from which planets are ejected due to fly-bys. 
/eject- The fractions listed were calculated for clusters with an 
initial half-mass radius of 0.38 pc, and only for encounters with 
mi = 0.6A/q . The calculated fractions assume that all solar-mass 
stars host a planetary system of the given type. In order to calcu- 
late the fraction of all stars affected, the values listed below must 
be multiplied by / p , the fraction of stars which hosts planetary 
systems. 
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interested in here, is to understand what fraction of plan- 
etary systems are significantly affected by fly-bys. In such 
systems, planet-planet scattering often leads to the ejection 
of one or more planets. Hence, the fraction of systems which 
eject planets is a good measure on the fraction of systems 
significantly changed by fly-bys. 

In the context of the extrasolar planet population, we 
are interested in the properties of the planets that are left 
orbiting stars after fly-bys. We describe some of these prop- 
erties in the second part of the section. 

9.1 Fraction of planetary systems damaged by 
fly-bys in stellar clusters 

We calculate the fraction of planetary systems from which 
planets are ejected, both immediately and in total, as well 
as the fraction of systems from which one or more planets 
is captured by the intruder in fly-bys as follows. First, we 
measure the number of stars in our stellar cluster simula- 
tions that have undergone fly-bys with r m i n < 100 au, TVfb. 
We have included the effects of multiple fly-bys in the cal- 
culation in a simplified way, by assuming that a star which 
has had iVfb fly-bys, is equivalent to iVfb stars having had 
one fly-by each. By doing this we neglect the combined ef- 
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feet of multiple fly-bys on one planetary system. Once we 
know iVfb, we can multiply this number by the fraction of 
fly-bys having r m i n < 100 au that lead to ionisation, cap- 
ture and /or the ejection of at least one planet in 10 8 years 
after the fly-by. Here, we only include fly-bys with intruder 
stars of mass mi = O.GMq. Having obtained the number 
of planetary systems, we then divide this by the number 
of stars hosting planetary systems in the cluster, to obtain 
the fraction of planetary systems which involved in ionising 
encounters, /ion, the fraction of system involved in fly-bys 
which lead to the capture of at least one planet, / cap , as well 
as the fraction of planetary systems that eject planets in 10 8 
years after fly-bys, /eject- 
In Table [5] we list /j on , /cap and /eject in clusters with 
r Ki = 0.38 pc for the 4G, GM, 4J and RlOa systems. 

The values of /i on , / cap and / c j cc t listed in Tableware 
lower limits. When calculating them we used the fraction of 
fly-bys causing ejections and/or capture of planets as cal- 
culated for mi = 0.6Mq. However, as can be seen in, for 
example, Table 1, fly-bys with higher mass stellar intruders 
are more damaging to planetary systems. We can, for the 4G 
system, include the increased damage done by fly-bys with 
higher mass as follows. We divide the encounters into three 
mass-bins (O.2M < mi < 0.8M Q , O.8Af < mi < 1.2M Q 
andl.2M0 < mi < 5Mq) and calculate what fraction of 
solar-mass stars have had an encounter within r m j n < 100 
au for each mass bin. We then multiply this number with 
the fraction of such fly-bys which cause the ejection of at 
least one planet in 10 8 years after the fly-by. We find that 
/eject for the 4G is then increased, from 0.11 to 0.14 n our 
reference cluster. 

In Fig. 13 we plot the fraction of systems, /eject, from 
which planets are ejected in 10 s years after a fly-by, as a 
function of the number of stars, TV, in the clusters. We in- 
clude the standard deviation of /eject as calculated for the 
4G system. This uncertainty is the standard deviation from 
the mean of the fraction of stars which undergo encounters 
in our cluster simulations. 

The solid line in Fig. 13 is a least-squares fit to /eject, 
as measured for the 4G system, where we assume that: 



/ 



eject 



oc N~* 



(9) 



The best fit gives 7 = 0.52 ± 0.11. 

The fraction of stars, /one, which encounter another star 
within a distance r m i n , during the lifetime of a stellar cluster, 
r, can be estimated as ([Heggie fc Hudl2003l l: 



/enc 



novdt 



(10) 



where n is the number density, a is the cross-section of the 
encounter and v is the velocity dispersion. 

The velocity disp ersion in a cluster de pends on its ra- 
dius, r h and mass, M (|Heggie fc Hut|[2003l '): 



However, all the clusters in our simulations have roughly 
the same size during most of their lifetime (rh ~ 2 pc), and 
hence we can write: 



u oc TV 



(12) 



The cross-section, a goes as a oc v~ for a given value of 
r m in when dominated by gravitational focussing (see Equa- 
tion and hence it scales with TV as: 



a oc — oc — 

v 2 N 



(13) 



while the number density, n scales linearly with the number 
of stars, n oc TV, as rh is constant. 

The lifetimes of clusters, r, increase with the number 
of stars, TV. In all the clusters, 50 per cent of fly-bys occur 
in the first 10 7 years, and 95 per cent of fly-bys in the first 
10 8 years. As an order of magnitude estimate, we can thus 
write: 



/enc 



navdt 



(14) 



where to = 10 , as the lifetime of all the clusters we have 
simulated are longer than 10 8 years. Combining these equa- 
tions we thus find: 



/enc OC TV 



1/2 



(15) 



We assume that the fraction of fly-bys within a given 
value of r m i n , which damage planetary systems, does not 
change with TV. Then we can write: 



/eject OC TV 



1/2 



(16) 



This result agrees well with our TV-body simulations. 
That /eject depends only weakly on TV is particularly inter- 
esting in terms of the the cumulative probability, P(N), that 
a given star forms in a cluster of size TV. In the sola r neigh- 
bour hood, observational surveys suggest that (see lAdamsl 
120101 . and references therein): 



P(TV) oc In AT 



(17) 



for TV = 100 — 2000. This means that all of our clusters 
simulated will contribute roughly equally to the stellar pop- 
ulation. Given that /eject varies only weakly with TV, low- TV" 
clusters can thus play an important role in populating the 
solar neighbourhood with planetary systems which have un- 
dergone planet-planet scattering caused by fly-bys. 



9.2 Properties of the planet population 



2 M TV 

V oc — oc 

rh r h 



(11) 



where we in the last step use the relation M = m x TV , with 
m the average mass of stars in the cluster and TV the number 
of stars. 



In terms of understanding how fly-bys can affect a given 
population of planetary systems, we must look at the prop- 
erties of these planetary systems long after the fly-by. In this 
section we discuss the eccentricities and semi-major axes of 
the planetary systems left around the host star 10 8 years 
after a fly-by. 
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Figure 14. The eccentricity distribution of planets with semi- 
major axisin the range 0.1 < a < 4.5 au, both observed (detected 
using radial velocities only, long-dashed line) , and from our simu- 
lations of fly-bys involving the 4G (solid line), GM (short-dashed 
line) and 4 J (dotted line). 



many planetary systems must contain planets with semi- 
major axes less than 4.5 au initially (i.e. the planets must 
have migrated there in a proto-planetary disc). The ob- 
served population thus contain such unperturbed systems 
with planets on rather circular orbits while the 4G popula- 
tion in Fig. [14] does not. 

It is interesting to note that eccentricity distribution 
of planets in the post-fiy-by 4G systems (Fig. I14|) differs 
from those which we find when simulating these planetary 
systems in binaries, where hierarchical systems, like the 4G, 
produced m uch fewer eccentric system s than the observed 
distribution (|Malmberg fc Daviesl f2009). The cause of this 
is that in fly-bys, the eccentricity of the inner planet is not 
only excited by planet-planet scatterings, as is the case in 
binaries where the Kozai mechanism operate, but also by 
direct interaction with the intruder star in very tight fly- 
bys. 



9.2.2 Semi-major axes 

In planet-planet scattering, the total energy of the plan- 
etary system is conserved. We can estimate the smallest 
semi-major axis which a planet can have after planet-planet 
scattering. The total binding energy in a planetary system, 
£"bind, depends only on the mass of the planets, m p ,i, their 
semi-major axes, a Pi i and the mass of the host star, Mh: 



9.2.1 Eccentricities 

In Fig. [T?] we plot the eccentricity distribution of the ob- 
served planets detected using radial velocities (long-dashed 
line). We also plot the eccentricity distribution of 4G (solid 
line), GM (long-dashed line) and 4J (dotted line) plane- 
tary systems as well as the system of run 10, series A from 
|Levison et al.|[l998r ) (dash-dotted line). The planetary sys- 
tems are "observed" 10 8 years after the fly-by, i.e. at the 
end of our simulations. When producing these distributions, 
we have included all planetary systems in our fly-by simu- 
lations (i.e. also those from which no planets were ejected 
in 10 s years). However, we only include planets on orbits 
with a semi-major axis less than 4.5 au, which is roughly 
the current detection limit in radial velocity surveys. Fur- 
thermore, we have limited the sample to planets on orbits 
wider than 0.1 au. Planets on tighter orbits than this will 
be tidally circularised, leading to a pile up of planets with 
e = in the observed sample, which is not representative 
of the mechanism which excite the eccentricities. All of the 
planetary systems which we have simulated here, have gas 
giants on orbits with semi-major axis larger than 4.5 au if 
no scattering has occurred. 

As can be seen, the eccentricity distribution of the 4G 
post-scattering systems is rather similar to the observed 
eccentricity distribution. Using the Kolmogorov-Smirnov 
statistic, which determines how different two sets of data 
are, we find that the two samples are consistent with origi- 
nating from the same underlying population (Pks = 0.61). 

We caution the reader that while the 4G eccentricity 
distribution appears to fit the observed planets rather well, 
it is important to note that we have excluded any contri- 
bution from unperturbed planetary systems in our data. In 
order to explain the observed semi-major axis distribution, 



-Ebind = — 



GA/ H m p ,i 
2a D . ; 



(18) 



where n p i is the number of planets in the system. 

Assuming that all but one planet is ejected in planet- 
planet scatterings, the semi-major axis of the remaining 
planet, at, is: 



ai = 



?7lf 



Ei=l m P,i/ ffl P,i 



(19) 



where itif is the mass of the remaining planet. Here we have 
assumed that the kinetic energy of the ejected planets is 
negligible. 

In, for example, the 4G system, where most often 
Jupiter is the remaining planet after scattering, we find 
a,f = 4.4 au, assuming that the planets were on their original 
orbits at the onset of scattering (i.e. the fly-by only changed 
the eccentricities of the planets). In more democratic sys- 
tems, the semi-major axes can be decreased by a larger fac- 
tor. In, for example, 4J equation ()19[) gives at = 2.6 au. 
However, many extrasolar planets are found on orbits much 
tighter than 2.6 au. Hence, although the eccentricities of 
the post-fly-by systems in Fig. [14] match the observed dis- 
tribution, their semi-major axes distribution does not. The 
observed planets on tighter orbits can, however, be produced 
by planet-planet scattering if the semi-major axes of the gas 
giants are significantly smaller than in the solar system at 
the onset of planet-planet scattering. We are therefore lead 
to conclude that migration must be working to bring gas 
giants in to distances of about 0.5 - 2 au before the sys- 
tem becomes unstable on its own or via the perturbations 
of passing stars or stellar binary companions. 
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10 DISCUSSION 

10.1 Fly-bys involving binary intruders 

When calculating the rate at which planets are ejected from 
planetary systems one should remember that the intruder 
may be a binary system. Such encounters may be more dam- 
aging to planetary systems th an fly-bys of single stars (see 
also I Adams fe Laughlin| [2001). Currently, we have included 
such encounters by setting the mass of the intruder equal to 
the sum of the masses of the binary components. However, 
by doing this, we neglect the possibility that these encoun- 
ters may lead to the formation of a transient stellar 3-body 
system, which could potentially be very damaging to any 
planetary system around the single star. 

From our JV-body simulations, we find that in about 1 
in 15 fly-bys which involve solar-mass stars, the intruding 
star is a binary system. As one third of stars in our clus- 
ter simulations are in binaries, one might expect the rate 
to be significantly higher. However, binaries are destroyed 
in encounters between two binaries, thus the po pulation of 
binar ies is depleted over time (see, for example, iHut et ail 
1 19921 ). 

Given the the binary population is really rather small in 
our clusters for most of the time, the contribution made to 
the fly-by encounters is small. If further, we assume that all 
encounters involving a binary lead to the ejection of at least 
one planet then we find that /eject for the 4G, as measured 
at 10 8 years, increases from 0.11 to 0.12. 

However, if the primordial binary fraction in the clus- 
ters were higher, which is suggested by observations (e.g. 
iDuauennov fe Mavoj[l99lh . the rate of fly-bys involving bi- 
naries would also be much larger. This in turn would cause 
a larger fraction of planetary systems to be strongly affected 
by fly-bys in clusters. We will present the results of simula- 
tions of clusters with a higher primordial binary fraction in 
a forthcoming paper. 

10.2 The effects of fly-bys on protoplanetary discs 

In this paper we only consider the effects of fly-bys after 
the phase of giant-planet formation is complete. However, 
an important aspect of fly-bys is also how they affect the 
planet-formation process. 

A fly- by o n a protoplanetary disc c an, for example, trun- 
cate the disc. iKobavashi fc Ida (|200lj ') find, using both an- 
alytical calculations and numerical simulations, that discs 
can be truncated at a radius rd = r m i n /3 in encounters 
with other stars, thereby halting planet formation at radius 
larger than rd. Hence, if, for example, a star has a fly- by 
with r m i n w 100 au, the disc will be truncated at a radius of 
about 30 au, the semi-major axis of Neptune. 

Fly-bys can also affect the properties of planets which 
form through core-accretion. Simulations of parabolic fly- 
bys of discs with on-going planet formation, suggest that 
fly-bys may cause planets to migrate outwards, and in crease 
the final mass of the planets |Fragner fc Nels on 2009). 

It is also possible that fly-bys can trigger instabilities 
in protoplanetary discs, inducing pla net formation via gray - 
itational instability. Simulations by iForgan fc Rice] (|2009l ) 
suggest that this is, however, not the case and that fly-bys 
instead tend to halt any fragmentation in the disc, thereby 
decreasing the rate at which planets form. 



The high number density of stars in stellar clusters 
means that radiation from nearby stars can affect the planet 
formation process via photo-evaporation. Photo-evaporation 
of discs decreases, or halts, planet formation. The rate at 
which this o c curs i n emb edded clusters was studied by 
lAdams et all (|2004 120061 '). They find that in 10 million 
years, photo-evaporation in a cluster containing about 1000 
stars can truncate protoplanetary discs down to a radius, rd, 
of rd = 36 au (Mh/Mq), where Mh is the mass of the planet 
host star. Thus, in such a cluster, photo-evaporation could 
affect giant planet formation, in particular around stars of 
somewhat lower mass than the Sun. 



10.3 Planetary systems in primordial binaries 

In our simulations of stellar clusters l|Malmberg et al]|2007T l , 
one t hird of stars are in primor dial binaries. Observations 
(e.g. IDuauennov fc Mavoij[l99ll ) suggest that the actual bi- 
nary fraction may be much higher than this, as about 50 per 
cent of stars in the field are observed to be part of a binary 
system. Increasing the binary fraction in our simulations 
would increase the rate of exchange encounters, thereby in- 
creasing the fraction of planetary systems which spent time 
in a binary system. 

In this paper we only consider the effects of encoun- 
ters on planetary systems around initially single stars. 
It is however quite possible that pla netary systems also 
form in primordial b i nary s ystems (e.g B arbieri et alj|2002l ; 
iThebault et~aH |2004 feood ) ■ While the inclination between 
the binary companion and the planets in such a system may 
initially not be large enough for the Kozai mechanism to 
operate, fly-bys can change t he or ientation of the compan- 
ion star. IParker fc Goodwin! (|2009T ) simulated the evolution 
of young embedded clusters for 10 Myr, with all stars in 
primordial binaries (i.e. the binary fraction equalled one). 
They found that 20 per cent of binary systems (which were 
initially all assumed to be co-planar) attained an inclination 
high enough for the Kozai Mechanism to operate. This frac- 
tion may become even higher in long-lived clusters, such as 
those studied here. We will re-visit this issue in a forthcom- 
ing paper. 



10.4 Fly-bys involving planetary systems with 
planets on tight orbits 

As mentioned in Section 9.2, fly-bys involving planetary sys- 
tems with similar properties to the solar system (i.e. the sys- 
tems discussed in this paper), cannot reproduce the semi- 
major axes distribution of the observed planets, / D b s (a). 
To explain the shape of / bs(<j), the primordial population 
of planetary systems (i.e. the population of planets before 
planet- planet scattering), must contain systems with planet 
on orbits significantly inside the ice line. Hence, many plan- 
ets must have undergone disc migration. 

In the current study we have not simulated the evolu- 
tion of such tighter systems involved in stellar fly-bys. As 
such, we are not currently able to quantify how the rate at 
which such systems are significantly changed in stellar clus- 
ters would differ with respect to solar-system-like planetary 
systems. However, we know that two competing effects will 
play an important role: 
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(i) Increasing a oute r /r m i n : For a given fly- by, the ratio be- 
tween the semi-major axis of the outermost planet, a ou ter, 
and r m in determines the magnitude of the perturbation on 
the planetary orbits. As such, one might expect that the in- 
stant effect of fly-bys involving tight planetary systems will 
be smaller than for fly-bys involving solar-system-like plan- 
etary systems. 

(ii) The stability of tight planetary systems: A system 
which has undergone disc migration may be more tightly 
packed than a solar-system-like planetary systems. As such, 
the perturbation needed in such a system to trigger planet- 
planet scattering might be considerably smaller. Their sta- 
bility may also crucially depend on two or more planets be- 
ing in a mean-motion resonanc e with each other (see, for 
example. iBarnes fc Q uinn 2004|). Even a rather small per- 
turbation, which breaks the mean-motion resonance, might 
then be enough to trigger planet-planet scattering. 

Hence, it is not clear how the effects of fly-bys, and 
how they affect the properties of the primordial population 
of planets, will change if many planetary systems contain 
planets which have migrated onto much tighter orbits. We 
will study this further in a forthcoming paper. 



10.5 The formation of extrasolar planets 

In Fig. fT5] we show an overview of the possible evolution- 
ary scenarios of planetary systems. The boxes show differ- 
ent types of systems, and the arrows identify how these may 
evolve. Planetary systems resembling the solar system, like 
those simulated in this paper, belong to the category "Plan- 
ets on wider orbits" . 

Black arrows identify an evolutionary scenario occur- 
ring without influence from any external bodies, while grey 
arrows identify evolutionary scenarios triggered by pertur- 
bations from other stars. 

For example, a "Protostar with disc" can evolve into 
a system resembling the solar system, either through nor- 
mal evolution (e.g core-accretion) or through fly-by induced 
fragmentation of the disc (e.g. gravitational instability). 

Such systems may remain stable, like the solar system, 
or dynamically evolve. If the system is long-term stable, as 
is the case for four gas giants, a fly-by or exchange encounter 
can turn it into an "Unstable system" . If the planetary sys- 
tem is initially not stable, as is the case of the four Jupiters, 
the systems will pass into the box labelled unstable systems. 

Unstable systems undergo planet-planet scattering, pro- 
ducing the types of extrasolar planetary systems which we 
observe today in radial velocity, transit and imaging surveys. 

One interesting possibility is that planetary systems 
which undergo exchange encounters will sometimes find 
themselves in extremely inclined binaries. If so, the Kozai 
mechanism can excite the eccentricities of planets to values 
close to unity. Planets will then have a very small perias- 
tron distances, causing them to tidally interact with the 
host star. The tidal interaction causes the planets' orbits 
to be circularised, greatly decreasing their semi-major axes. 
This behaviour, often called Kozai migration, can lead to 
the formation of hot Jupiters f rom planets on initially muc h 
wider orbits (see, for example. iFabrvckv fc TremainelfeoO'it) . 
A key point here is that Kozai migration will often leave the 
hot Jupiters significantly inclined with respect to the rota- 



tional axes of the host stars. For transiting hot Jupiters, it is 
possible to measure the projected inclinations of the planets' 
orbits with respect to the stars' rotati onal axes, the so - called 
Rossiter-McLaughlin effect. Recently, iTriaud fc et~ai1 (|2010h 
showed that the distribution of inclinations for a large set 
of observed hot Jupiters follows that expected if they were 
all made through Kozai migration, and is inconsistent with 
formation via e.g. type II disc migration. 

Planet-planet scatterings will ordinarily place planets 
on transient wide orbits (i.e. with semi-major axes greater 
than 100 au) in both the intrinsically unstable systems and 
in those systems where instabilities are induced by encoun- 
ters. Planetary systems which undergo scatterings are in 
general only a few 100 million years old. Hence, if planet- 
planet scattering is responsible for the eccentricities that we 
measure in the observed exoplanet systems, planets on wide 
orbits will be o bserved in imaging surveys of nearby solar- 
type stars fe.g. lHeinze et al]l2010h . 



10.6 Observational signatures of fly-bys 

As shown in this paper, the principal effect of fly-bys is 
that they may trigger planet-planet scattering in other- 
wise stable planetary systems. However, the post-scattering 
systems cannot easily be distinguished from systems in 
which the scattering was trigg ered because the systems 
were intrinsically unstable (e.g. iMoorhead fc Adams! [20081 ; 
lJuric fc Tremain j|200d IChatteriee et al.ll200Sft ~ 

One important aspect of fly-bys is, however, that only a 
fraction of stars suffer fly-bys. Hence, a large fraction of stars 
are singletons, and any planetary systems around them have 
not been perturbed by other stars. Hence, we predict that 
once planets can be detected in the semi-major axes region 
between 5 and 30 au around other stars many systems with 
similar properties to the solar system should be observed, if 
most systems start out with similar properties to the solar 
system. 

Furthermore, as stellar encounters occur throughout the 
entire life-times of stellar clusters, we would expect that the 
fraction of planetary systems which are perturbed increase 
with cluster age. Hence, the properties of the population 
of planets in stellar clusters will depend on the age of the 
cluster. 

Another important observational signature of planet- 
planet scattering, whether it is triggered by interactions with 
ot her stars or not, i s , as is described both in th is paper and 
in IScharf fc Menoul (120091 ); IVeras et al.l i|2009h , that it will 
produce planets on very wide orbits. Such planets must be 
observed in future imaging surveys, if planet-planet scat- 
tering is responsible for the observed eccentricities among 
extrasolar planets. 



11 SUMMARY 

We have numerically simulated fly-bys involving planetary 
systems resembling our own solar system. Our motivation 
was to understand how a population of planetary systems 
evolve in stellar clusters, the birth environment of a large 
fraction of all stars. We have combined the fly-by simulations 
with our previous simulations of young stellar clusters, to 
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Figure 15. Overview of the evolution of planetary systems. Grey arrows identify interactions with other stars, most of which occur in 
young stellar clusters. Examples of such are fly-bys and exchange encounters involving binary systems. Planetary systems similar to our 
own solar system belong to the box labelled "Planets on wider orbits" . 



calculate the effect on the population of planetary systems 
in the field today. The main results of the paper are: 

(i) Close fly-bys can lead to the immediate ejection of 
one or more planets from the host system. Such planets 
are either left unbound (ionisation) , or are captured by the 
intruder star (exchange). Captured planets are often left 
on wide and moderately eccentric orbits, leaving them de- 
tectable in imaging surveys. Planets which remain bound to 
the host star after a close fly-by are most often on signif- 
icantly more eccentric orbits. In a moderately-sized stellar 
cluster (TV = 700 stars, — 0.38 pc) about 3 per cent of 
solar-system-like planetary systems will be involved in fly- 
bys leading to ionisation and/or capture of planets. 

(ii) Wider fly-bys, which do not immediately remove plan- 
ets from the host star, may still cause the ejection of planets 
on longer time scales by planet-planet scattering. The time 
between the fly-by and the onset of planet-planet scattering 
can be very long. For the four gas giants of the solar system 
we find that we can accurately predict if the system will 
suffer planet-planet scattering within 10 8 years after a fly- 
by (e.g. Table 1), while the same is not easily done for the 
system where the masses of the planets have all been set to 
that of the planet Jupiter (4J) (e.g. Table 3). Planet- planet 
scattering significantly increases the fraction of planetary 



systems from which planets are ejected due to fly-bys. In, 
for example, the case of the 4G planetary systems and en- 
counters with r m i n < 100 au and mi = 0.6Mq the fraction 
of systems which suffer ejections is increased by a factor of 
three from the time of the fly-by until 10 8 years later. 

(iii) If the intruder in a fly-by is a very low-mass star 
or a brown dwarf, a sufficiently close fly-by may leave the 
intruder bound to the host star, potentially leading to dras- 
tic changes in the orbital properties of the planets orbit- 
ing the star. This can significantly increase the effect that 
such fly-bys have on planetary systems compared to what 
we would naively expect; for stellar mass intruders the mag- 
nitude of the change in the planetary orbits due to the fly- 
by decreases with decreasing intruder-star mass, suggesting 
that sub-stellar mass intruders would do very little damage 
to planetary systems. We have performed a large number of 
scattering experiments between a single-planet system and 
low-mass intruders to measure the rate at which the intruder 
may become bound to the host star. We find that, depending 
on the velocity dispersion of the clus ter, and given the log- 
normal IMF of lChabrierl i|200ll . l2002r ), up to few per cent of 
solar-system-like planetary systems may be involved in fly- 
bys that result in the capture of the low-mass intruder. Of 
course, if real clusters have an e ven larger number of low- 
mass intruders than suggested bv lChabrierl |200ll . [2"002l ). the 
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fraction of systems involved in encounters with such objects 
will increase. 

(iv) Fly-bys in general excite the eccentricities of plan- 
ets left bound to either the host or the intruder star, both 
through direct scattering off the intruder, and through any 
subsequent planet-planet scattering. We find that the ec- 
centricity distribution of planets which have been scattered 
inwards due to fly-bys in the 4G system is similar to that of 
the observed extrasolar planets (see Fig. 7 and discussion in 
section 9.2). Planetary systems where the masses are more 
similar (e.g. GM & 4J) tend to overproduce eccentric planets 
compared to the observed extrasolar planets. 

(v) Planet-planet scattering, triggered either by fly-bys 
or by primordial instability decreases the semi-major axes 
of one or more planets while leaving others on either much 
wider orbits, or unbinds them entirely. In our simulations of 
the four gas giants, we find that planets are only rarely scat- 
tered onto orbits tighter than 4 au, which is not consistent 
with the semi-major axes of observed planets. To reproduce 
the extrasolar planet population it is thus needed to allow 
for disc migration before the fly-by, allowing for planets to 
be on orbits inside the ice-line. 

(vi) During the phase of planet-planet scattering, one or 
more planets are successively scattered onto wider and wider 
orbits until they become unbound. Such planets can, since 
they are so far from the host star, be observed using imaging 
techniques. 

While not studied in this pap er (but see 

iMalmberg et al.1 120071 : iMalmberg fe Daviesl l2009h another 
important type of encounters is binary-single encounters, 
in which a single star encounters a stellar binary system. 
Such encounters can, depending on the masses of the stars, 
leave the initially single star in the binary, replacing one of 
the original binary components. The key point here is that 
in such binaries, any planetary system will be randomly 
oriented with respect to the stellar companion, and as such 
will often be highly inclined with respect to the latter. In 
highly inclined binaries (i > 39.2°) the Kozai mechanism 
operates, exciting the eccentricities of the planets and 
thereby potentially triggering planet-planet scattering. 

We have combined our results of how fly-bys affect plan- 
etary systems with the rate at which encounters in young 
stellar clusters occur, to measure how important such en- 
counters are. As an example, we can look at the effects of 
fly-bys on the four gas giants of the solar system. We find 
that in between 5 and 15 per cent of such systems, fly-bys 
will trigger planet-planet scattering, leading to the ejection 
of one or more planets. The fraction of systems from which 
planets are ejected varies slowly with the number of stars 
in the cluster, N. More stars are believed to form in low- AT 
clusters than in more rich clusters. Hence, the contribution 
from rather small clusters to the field population of plane- 
tary systems that have been made unstable by fly-bys, can 
be significant. An important point regarding the fraction of 
systems that are affected by fly-bys concerns the properties 
of the systems before any stellar encounters. To explain the 
observed semi-major axes distribution many systems must 
initially have contained gas giants on much tighter orbits 
than what is seen in the solar system. The effects of fly-bys 
on such systems have not yet been explored, but will be the 
topic of a forthcoming paper. 



In summary, fly-bys in young clusters will significantly 
change the properties of the population of planetary systems 
which initially resemble the solar system. 
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APPENDIX A: ANALYTICAL FORMULAE 
PREDICTING THE EFFECTS OF FLY-BY'S 

In order to write down the equations we use the standard 
Keplerian orbital elements describing the shape of the orbits, 
as well as the orientation of the intruder star's orbit with 
respect to the orbit of the planet. 

Let, a p be the semi-major axis and e p the eccentricity 
of the planet long before the encounter while ei is the ec- 
centricity of the intruder star in its hyperbolic orbit around 
the barycentre of the planetary system. 

In order to define the orientation angles we define the 
line of nodes as the line formed by the intersection of the 
orbital planes of the planet and the intruder star. The as- 
cending node is the point on this line at which the intruder 



star cross the orbital plane of the planet, measured in the 
orbital plane of the planet in the the same direction as the 
angular momentum vector of the planetary system. The lon- 
gitude of the ascending node, Q is measured from the peri- 
centre of the planet in its orbital plane. Furthermore, i is 
the inclination between the two orbital planes and ui is the 
longitude of pericentre of the intruder star, measured in its 
orbital plane from the ascending node, in the direction of its 
motion around the planetary system. 

Finally, we define tuh to be the mass of the host star, 
m p to be the mass of the planet, mi to be the mass of 
the intruder star and r m i n to be the minimum separation 
between the stars in the fly-by. We write the sum of the 
mass of the planetary system (planet and host star) as M ps 
and the sum of the mass of the planetary system and the 
intruder star as Mtot- 

For planets o n initially eccentric o rbits the change in 
eccentricity, Se is l|Heggie fc Rasio 1996|): 



8e = 



15 



mi 



y/M ps M tot \r 



3/2 



(l + e : )3/ 2 



x < sin i sin 2Q 



arccos(— 1/ei) + -y/e^" 



+ 



r 2 
— [(1 + cos i) cos 2cj sin 2Q- 



2 cos i sin 2lu cos 2fi] 



N3/2 



(Al) 



This result vanish if the eccentricity of the planet equals 
zero, which is why one must include higher order terms in the 
derivation. For a distant encounter the chan ge in eccentricity 
for a planet on an initially circular orbit is ( Heggie fc Rasiol 
Il996h : 



5e 



15 mi \nin — m p \ ( M p 



liJ -ps 



Mu 



1/2 / \ 5/2 

/ a p N 



e?(l + ei) 



5/2 



A(ei)(l- -j sin i) 

1 2 



+ — (e 2 — 1) J / 2 (1 — 5 sin 2 u) sin 2 i 



15 



+ COS LO 



/ 1 (ei)(l- i sin 2 i) 



i 2 2 -i \5/2/-, r • 2 ■ 2 -n 

+ — (ei — 1) ' (1 — 5sm ojsm i) 



1/2 



(A2) 



where 



/i(ei) = e\ arccos(— 1/ei) + 



15 



(-2 + 9ei + 8ef) (A3) 



Finally, for tight encounters involving planets on 
initially circular orb its, the change in eccentricity is 
jHeggie fc Rasiol 1 1996T ): 
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5e = 3^?^(^V /4 ^+ - 1 ) 374 



Ml] 



x exp 



1/2 



r m i n ^ 3 / 2 v^l ~ 1 _ arccos (1/ei) 
(ei - l) 3 / 2 



X cos — 



cos — I — sin — 
2 9 2 



+ — cos g 2 ° OS 



1/2 



(A4) 



In fly-bys involving planetary systems containing mul- 
tiple planets, not only the change in eccentricity, e, but also 
the change in inclination, i, is important. Thus, we have, us- 
ing the same approach as in lHeggie fc Rasiol (| 19961) . derived 
an equation for the change in inclination of the planet, Si p , 
with respect to the plane of its original orbit: 



. . . 37rmi 
sm [5i p ) = — — 



3/2 

ip • 



nun 



MtotMp S (l-e|) \r 
x sin i cos i [(1 + 3e 2 ) 2 sin 2 0, + (1 — e 2 ) 2 cos 2 Ql 1//2 (A5) 

This expression is for a fly-by with the intruder on a 
parabolic orbit and has been averaged over the argument of 
pericentre (u). 
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